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Abstract. - Atomic Force Microscope images of a crack intersecting the free surface of a glass 
specimen are taken at different stages of subcritical propagation. From the analysis of image pairs, 
it is shown that a novel Integrated Digital Image Correlation technique allows to measure stress 
intensity factors in a quantitative fashion. Image sizes as small as 200 nm can be exploited and the 
surface displacement fields do not show significant deviations from linear elastic solutions down 
to a 10 nm distance from the crack tip. Moreover, this analysis gives access to the out-of-plane 
displacement of the free surface at the crack tip. 
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1 . The mechanisms of subcritical crack propagation in sil- 
^^icate glasses have been the subject of extensive research 
(cf. [UH] for recent reviews). The development of ad- 
vanced AFM techniques to probe in- situ crack propagation 
f-^or post-mortem crack surface morphologies have recently 
OOled to explore these mechanisms at their relevant nano- 
r ~"J metric scale leading to remarkable observations on plastic 
CN crack tip damage [3]|4], stress-induced ion migration [5] 
and capillary condensation inside the crack tip cavity [6] . 
However, a proper understanding of these observations in 
• • terms of sound mechanical modeling has been hampered 
. ^ by the lack of a technique to measure the stress and strain 
^fields at the crack tip with a nanometric resolution, and 
^nthe proposed interpretations are still subject of debate [2]. 
5^ Another weakness of the in-situ AFM observations is their 
limitation to access the external surface of the sample, that 
is intersected by the propagating crack front with a small 
tilt angle, thus making the local three-axial condition of 
stress not trivial 0. The present letter reports a con- 
comitant solution of all these problems. By extending the 
Digital Image Correlation technique to AFM topograph- 
ical images of in-situ crack propagation in a silica glass 
it was possible (1) to characterise properly the 3D sur- 
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face displacement fields with nanometric resolution, (2) to 
show that they can be fitted by appropriate elastic solu- 
tion down to a 10 nm distance from the crack tip and (3) 
that the locally measured surface stress-intensity factor 
is in agreement with the macroscopically measured value 
within 5% uncertainty. 

Digital Image Correlation (DIC) is a technique which 
allows one to measure displacement fields by matching 
a reference with a deformed image, most often an op- 
tical image. Introduced a long time ago in the field of 
solid mechanics by Sutton et al. [81, this technique has 
known a very rapid development due to the wide avail- 
ability of digital imaging devices, and the increasing per- 
formance and reliability of analyses. In particular for 
cracks, since the pioneering work of Sutton [9], recent 
works [TQHT2] have shown that stress intensity factors can 
be estimated very accurately either directly through In- 
tegrated DIC (IDIC) which incorporates analytic crack 
fields [TO], or through a tailored post-processing of the 
displacement field [T m[TT| [T3]. However, most applications 
of DIC deal with optical images, and very few with Atomic 
Force Microscopy (AFM) images. Among those, the most 
advanced published works concern thin polycrystalline sil- 
icon films [T4UT8] . In particular, the measurement of a 
uniform strain field from AFM images at different scales 
down to scale 1x2 fim 2 images [T7] allowed to estimate the 
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elastic properties of polycristalline Silicon used for MEMS. 
More recently, stress intensity factors were evaluated for 
cracks propagating in such films [18]. AFM was used to 
locate the crack geometry. Stress Intensity Factors (SIFs) 
were then indirectly estimated through an elastic finite el- 
ement simulation of the test and specimen geometry. In a 
different spirit, a recent study by Xu et al. [19] investigated 
systematic bias such as drift and distorsion of large AFM 
images (6.4 x 6.4 /im 2 ) on unstrained samples in order to 
estimate (and reduce) strain measurement uncertainties. 
These two sets of studies have demonstrated the feasibil- 
ity of using DIC on AFM images, mostly for simple strain 
fields and large scanned regions. 

The present work is aimed at a direct determination 
of stress intensity factors through DIC at a much lower 
scale than previously considered, (from lxl jim 2 down 
to 200 x 200 nm 2 , i.e. about 100 times smaller than in 
previously reported studies [H]), without any recourse to 
numerical finite-element simulations. This analysis also 
provides information on inelastic processes at play in the 
crack tip neighborhood through unresolved discrepancies 
after accounting for a simple linear elastic-brittle behav- 
ior. Applying DIC requires to overcome a number of very 
significant challenges: As the size of the region of inter- 
est is reduced, the noise level of AFM images becomes 
very significant. Moreover, images consist in topographic 
measurements. In contrast with classical speckle patterns 
which are simply advected by the displacement field, the 
topography is affected by out-of-plane motion of the free 
surface. In the present study, we will see that the typ- 
ical peak-to-valley roughness of say 1 /im size images is 
of order 3 nm. The out-of-plane displacement difference 
due to a crack at the critical stress intensity factor will 
be estimated below as being of order 2 nm, and thus can- 
not be ignored. Thus, the very foundation of DIC, i.e. 
the conservation of optical flow, has to be revisited and 
generalized to account for this effect. 

In the present experiments, fractures were initiated and 
propagated on a Double Cleavage Drilled Compression 
(DCDC) set-up (cf. Fig. [IJ at a constant temperature of 
(23 ± 1)°C in a leak-proof chamber under an atmosphere 
composed of pure nitrogen and water vapor at relative 
humidity (23 ±1)%. The DCDC test set-up is particu- 
larly convenient for these studies due to its excellent sta- 
bility and several theoretical modeling were devoted to 
it [2U1I2T]. The parallelepipeds DCDC samples (4 x 4 x 40 
mm 3 , ± 10 /im) of pure fused silica glass (Suprasil 311, 
Heraeus, Germany) were polished to a RMS roughness 
of 0.5 nm (for an area of lxl /im 2 ) and a hole of radius 
R = (500 ± 10) /im was drilled at their center to trigger 
the initiation of the two symmetric fractures of length c. 
AFM observations are done in tapping mode on a D3100 
from Veeco Metrology Inc., Santa Barbara, CA. The de- 
tails of the setup and techniques can be found in [3j[6]. 
After an extensive stabilization consisting of several hours 
of stationary propagation and imaging, all the images pre- 
sented here were acquired within two-hours, so that crack 




Fig. 1: Experimental setup: (a) Sketch of the DCDC geometry; 
(b) picture of the experiment. 

propagation conditions can be considered as stable at a 
propagation velocity v = (0.7 ±0.1) nm/s and AFM drifts 
are minimized. For an applied force F = (1844 ±4) N and 
an average crack length c = (6145 ± 10) /im, the SIF can 
be estimated to Kj = (0.39 ± 0.02) MPa.m 1 / 2 according 
to Ref. [20]. 

Let us denote by f(x) and g(x) the topographic images 
of the reference and deformed states, where x are the co- 
ordinates in the observation plane. / and g are the height 
of the observed surface in nanometers. The generalization 
of the so-called "optical flow conservation" which relates 
the two images can be written 

g(x) = f(x + u(x)) + v(x) (1) 

where u (resp. v) denotes the in-plane (resp. out-of-plane) 
displacement field. 

In order to develop a robust method, we resort to an "In- 
tegrated" DIC (IDIC) approach as introduced in Ref. [TO] . 
The displacement field (u(x),v(x)) is searched for as a 
combination of a few expected fields $> n (x) 

10 

(u(x),v(x)) = a>i®i(x) (2) 

i=l 

Those fields consist in rigid body motions 1 < n < 6 and 
mode I crack displacement fields detailed in the following: 
The latter fields are well known in either plane stress or 
plane strain conditions [22], which unfortunately do not 
strictly apply for a crack front emerging on the free sur- 
face of a 3D sample [21]. Nevertheless, the in-plane com- 
ponent u(x) can be well approximated by the plane stress 
condition. We introduce the notation (f(x) for such an 
in-plane displacement field for a crack tip located at the 
origin, and crack path along the x < semi- axis. The 
displacement field normal to the free surface is not known 
analytically. The latter, for a unit mode I stress intensity 
factor, is denoted ip(x) for the same reference geometry. 
Because of their different status, the in-plane and out-of- 
plane displacements will be treated in the sequel as two 
independent fields. 

An additional difficulty is that as we follow a crack prop- 
agation under a constant load, the crack tip is present in 
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both images although at different locations. Let us call 
Xo and X\ the crack tip positions in / and g respectively. 
The displacement field to be applied to the reference image 
thus consists in closing the crack at point Xq and opening 
it at X\. Thus the 7 th and 9 th basis fields are cp(x — Xo) 
and (p(x — X\) respectively, while the 8 th and 10 th fields 
are the out-of-plane components ip(x — Xo) and ip(x — X\) 
respectively. For a steady state propagation under the 
same stress intensity factor, we can write a 7 = —ag and 
as = — cliq. Those equalities can be enforced using La- 
grange multipliers. Thus the problem consists in evalu- 
ating the above 10 unknowns, a n , (8 degrees of freedom) 
through a weak form of the extended optical flow conser- 
vation to be minimized: 

T= J (g(x) -f(x + ai ^\x)) - a^f\x)) 2 ^ 

(3) 

where <E>-^ and designate respectively the in-plane 
and out-of-plane displacement. When suitably normalized 
(using the pixel size and the elastic properties of the silica 
glass specimen), we see that the SIF is directly estimated 
as the a 7 parameter value. The problem as it is written is 
strongly non-linear because of the occurrence of unknowns 
as arguments of the topography /. To avoid secondary 
minima trapping the determination of the displacement 
is performed hierarchically from the minimization of the 
above T functional from coarse-grained images down to 
the original images. Each minimization is performed it- 
eratively resorting to a linearization in the unknowns a n . 
This procedure is similar to the one introduced by Hild et 
al. [1H. 

Finally, the above computation requires the knowledge 
of the crack tip location. If, at a large scale, the crack path 
is clearly visible on the topographic images, the exact po- 
sition of the tip is more uncertain (cf. Fig.[2fa)). However, 
the hydrophilic nature of the glass surface induces a capil- 
lary condensation of water at the crack tip, and when the 
AFM tip probes the immediate vicinity of the crack, the 
formation (and disruption) of a water meniscus induces a 
significant change in the phase angle measurements [6] as 
can be seen in Fig. [2^b). Thus, albeit phase angle maps 
are much more homogeneous than topographic ones, and 
hence can hardly be used for IDIC, the crack path and tip 
can be resolved precisely on those images. We should note 
that the present phase images are acquired using the repul- 
sive tapping mode, which allows for a resolution of 10 nm 
on the lateral positioning of the crack tip, but good results 
can also be obtained using either the attractive mode or 
a spatially mixed mode between the surface and the con- 
densate (cf. [IOCS])- In addition, starting from these obser- 
vations of the crack tips on both reference and deformed 
images, a final optimization can be performed based on 
the minimization of the functional T. 

One major difficulty is that the analytical form of the 
out-of-plane field ^{x) is unknown. The general problem 
of the emergence of a crack on a free face, is a difficult 
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Fig. 2: Topographic (a) and phase (b) images. The whole crack 
path appears to be clearly visible on the phase map, but the 
crack tip can hardly be perceived on the topographic image. 

issue. In the framework of linear elasticity, at the very 
point of intersection between the crack front and the free 
surface, i.e. the observed "crack tip", a point singularity 
occurs, whose exponent depends continuously on the angle 
of the crack front incidence to the external surface and on 
the Poisson's ratio [25 . A detailed theoretical and numer- 
ical study of this effect applied to glass and specifically in 
the very same DCDC test as the one here considered can 
be found in Ref. [21]. 

It has been argued [7] that the crack front geometry 
adjust itself at propagation so that the singularity expo- 
nent is the same, 1/2, as for a planar crack. Thus, as a 
rough approximation, we first used tjjo(x) = vT^l- ^1- 
though, not equal to the actual displacement field, yet 
it allowed to register both images, with the appropriate 
in-plane displacement. Once this is achieved, we can ex- 
tract the height difference Sz(x) between the reference 
and deformed surfaces, ignoring the ipo(x) contribution. 
This difference (or "residual") should be expressed as 
i(j(x — xo) — ^(x — xi) . A simple algebraic expression for ijj 
is proposed as a systematic Fourier expansion respecting 
the mode I crack symmetries 

ip(x) = y/\x\ a n cos(n<9/2) (4) 

n 

where is the polar angle of the current point x with 
respect to the crack tip. Low orders for the parameter < 
n < 2 are chosen. This allows for a direct experimental 
determination of the out-of-plane displacement field which 
is then used in the IDIC analysis. 

A series of 5 images of the same lxl /im 2 zone swept by 
a crack is first analyzed. An example of the residual map 
without out-of-plane displacement is shown in Fig. [3^a). 
From the set of the four residual fields, a least squares 
regression is implemented in order to identify the out- 
of-plane displacement based on the algebraic expression 
Eq. H a = (-0.39,-0.94,1.0) is obtained. Figure E^c) 
shows that no obvious long range trend is left in the final 
residual. The out-of-plane displacement due to the open- 
ing of a crack at the origin is shown in the same Figure. 
It is to be emphasized that the present analysis cannot 
determine the absolute reference plane and thus ip is de- 
termined up to a linear function in x. Yet, the order of 
magnitude of the displacement amplitude is quite compa- 
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Fig. 3: (a) Raw residual field obtained after correction of in- 
plane displacement. Note that the crack path is masked for 
the IDIC analysis; (b) Fitted displacement field, i/j(x — Xo) — 
ip(x — jci). (c) Remaining residual after fit subtraction. The 
mask has not been included to make the crack visible; (d) Out- 
of-plane displacement ip(x) associated with the crack opening. 
A cut through the symmetry plane y = is shown, with the 
crack tip located at the origin. 



rable to the result of numerical simulations reported by 
Fett et al. [21]. 

This expression of the out-of-plane displacement is now 
used in the sequel. The quality of the IDIC measurement 
is evaluated by measuring the normalized residual, 77, i.e. 
the standard deviation of the remaining residual, normal- 
ized by the peak-to-valley roughness of the original to- 
pographic image. In the present cases, even though the 
kinematic field was adjusted with very few parameters, 
the normalized residual 77 was very stable at about 1.6%, a 
remarkably low value which gives confidence in the conver- 
gence of the procedure. The measured SIF values (assum- 
ing a Young's modulus of 72 GPa, and a Poisson's ratio 
of 0.17) were successively Kj = 0.38, 0.36, 0.44 and 0.38 
MPa.m 1 / 2 . The stability of the four independent measure- 
ments gives an a posteriori indication on the uncertainty of 
the toughness measurement, 0.39 ±0.04 MPa.m 1 / 2 , which 
is in excellent agreement with the macroscopic value of 
0.39 ± 0.02 MPa.m 1 / 2 estimated independently. Note the 
huge difference of length scales used in both estimates 
(millimeters vs. nanometers respectively for the macro- 
scopic evaluation and present analysis). 

Going to the smaller scan size of 200 nm, a number of 
additional difficulties arises: First the precise reposition of 
the same zone is done within a few tens of nanometers, i.e. 
a large fraction of the scanned image and such a large off- 
set is more difficult to handle for DIC. This difficulty was 
taken care of by a prior rough determination of the mean 
translation. Second, the fluctuations in the "slow scan" di- 
rection are much more intense than in the fast scan direc- 
tion (here transverse to the cracks). This high frequency 
noise of amplitude 0.9 A is only partly corrected by the 




Fig. 4: Analysis of 200 nm size images, (a): Reference image; 
(b): Deformed image; (c): Residual map (where the masked 
region around the crack can be seen); (d): U x displacement 
map (in nm); (e): U y displacement map (in nm). (f): U z 
displacement map (in nm); 



"flattening" filter used and it makes the computation of 
image gradient insecure. To further limit this noise, a low- 
pass Fourier filtering was applied on the image at the cut- 
off frequency for which the slow and fast scan power spec- 
tra depart from each other. For 200 nm images, although 
being composed of 512x512 pixels, wavelengths below 4 
pixels (i.e. about 1.5 nm) were filtered out. Another ma- 
jor difficulty comes from uneven step sizes in the slow scan 
direction. A simple bare eye inspection does reveal that in- 
deed, images appear to be either dilated or expanded along 
the slow scan direction by a significant amount. Similarly, 
due to inaccuracies in the reposition at the beginning of 
each new line an artificial simple shear can be introduced. 
Thus, we implemented additional degrees of freedom in 
the kinematics, (i.e. uniform and linearly varying strain 
along the slow scan direction, and uniform shear strain) 
corresponding to these artefact deformation modes. Fi- 
nally, it is worth mentioning that pointing at the crack 
tips, even from phase images, is quite delicate and does 
contribute to the final uncertainty. 

With those additional procedures, image series could 
be analyzed successfully, in the sense that the residual, 
77, could be brought down to 5-6%. A visual check based 
on the comparison between reference and corrected de- 
formed images indeed confirms the reliability of the image 
matching, and hence the robustness of the algorithm, in 
spite of the significant level of noise in the images. Mea- 
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surements of the mean strain along the slow-scan direc- 
tion revealed strains which could be as large as 8%. In 
the example shown in Fig. HJ the SIF is estimated to be 
0.44 MPa.m 1 / 2 , the axis strain is 1.6%, and the resid- 
ual amounts to r] = 5.4%. Ignoring the most strained 
image pair, the SIF value was estimated in the range 
0.37 < Ki < 0.46 MPa.m 1 / 2 . Varying the size of ana- 
lyzed domain, and the cut-off frequency induces only very 
slight modifications in the estimates of the SIF. The esti- 
mate Kj = 0.41 ± 0.05 MPa.m 1 / 2 is again consistent with 
the macroscopic determination and the larger scale AFM 
IDIC estimate. Most of the uncertainty lies in the precise 
identification of the crack tip location. As can be seen in 
the residual map shown in Fig. H^c), no systematic bias 
can be distinguished even in the vicinity of the crack tips. 

Based on those quantitative observations, no evidence 
of a breakdown of linear fracture mechanics — eventu- 
ally due to micro-damage or plasticity — could be ob- 
served even in the immediate vicinity (within 10 nm) of 
the crack tip. The present IDIC technique has thus proven 
to be suited to measuring stress-intensity factors in brit- 
tle materials by comparing pairs of AFM images of the 
external sample surface taken at different stages of crack 
propagation. Moreover, this analysis gives an original in- 
sight into the out-of-plane displacement fields of the free 
surface in the neighborhood of a propagating crack tip, 
which is presently an open subject of research. 
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